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Predicting the Aoki Phase using the Chiral Lagrangian* 
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This work is concerned with the phase diagram of Wilson fermions in the mass and coupling constant plane 
for two-flavor (unquenched) QCD. We show that as the continuum limit is approached, one can study the lattice 
theory using the continuum chiral Lagrangian, supplemented by additional terms proportional to powers of the 
lattice spacing. We find two possible phase structures at non-zero lattice spacing: (1) There is an Aoki phase of 
spontaneously broken flavor and parity, with two massless Goldstone-pions, and a width Amo ~ a 3 4 ; (2) There is 
no spontaneous symmetry breaking, and all three pions have equal mass of order a. Present numerical simulations 
suggest that the former option is realized. 


Some time ago, Aoki proposed that lattice 
QCD with Wilson fermions has a phase in 
which flavor and parity are spontaneously bro¬ 
ken, thereby providing a dynamical explanation 
for the presence of massless pions at non-zero lat¬ 
tice spacingjl]. His proposed phase diagram is 
shown in Fig. [I]. This suggestion has been stud¬ 
ied using analytical and numerical methods, and 
somewhat contradictory conclusions have been 
drawnjj. We have recently shown that one can 
analyze the phase structure analytically close to 
the continuum limit Q. In this paper we summa¬ 
rize the salient points of our workQ 

We consider lattice QCD with two degenerate 
Wilson fermions. The action is invariant under 
flavor SU (2) and parity, but explicitly breaks all 
axial symmetries for any value of the bare mass. 
Aoki’s proposal is that the condensate (Ats 0 ^) 
is non-zero in region B (with cr 3 acting in flavor 
space), thus breaking the flavor symmetry to a 
U{ 1) subgroup and producing two massless Gold- 
stone pions throughout B. 

To investigate this suggestion analytically, we 
first write down an effective continuum descrip¬ 
tion of the lattice theory valid close to the con¬ 
tinuum limit, and then examine the long distance 
behavior of the effective theory using standard 
chiral Lagrangian techniques. The effective con¬ 
tinuum theory can be described by a Lagrangian 
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3 For detailed references to earlier work see Ref. [2J 

4 M. Creutz had previously performed a similar analysis 
based on the linear sigma model fe|. 
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Figure 1. The phase diagram proposed by Aoki: 
g is the gauge coupling and mo the dimensionless 
bare quark-mass. The continuum-like unbroken 
phases are labeled A, and the flavor and parity 
broken phase B. The continuum limit of partic¬ 
ular interest is that at mo = 0, g = 0. 

in which the usual terms have been supplemented 
by higher dimensional operators proportional to 
powers of the lattice spacing |J]. These additional 
pieces are constrained only by the symmetries of 
the lattice action, which means that they need 
not respect chiral symmetry, and their enumera¬ 
tion is identical to that carried out as part of the 
improvement program js|]. As will become appar¬ 
ent shortly, we are interested in physical masses 
of order the lattice spacing, i.e. m ~ aA * 2 , with A 
an abbreviation for Aqcd- In this mass regime, 
with a suitable choice of quark variables, the ef- 






fective continuum description of the lattice theory 
can be written^) 

£ e ff = £ g + V>QZ) + m)ip + biatpia^F^il) , (1) 

where £ g is the gluon Lagrangian, and we have 
dropped terms of 0(a 2 ). The dimensionless pa¬ 
rameter bi ~ 1 is a function of the coupling g, 
while the physical mass is m ~ (mo — m c (g))/a , 
with mo the (dimensionless) bare mass and rh c a 
critical massfi Both the mass term and the Pauli 
term in (pj) break chiral symmetry. A key point 
in the following is that when we adjust the bare 
mass at fixed g (and thus fixed a) we change the 
relative importance of these two terms. 

The next step is to write a generalization of the 
continuum chiral Lagrangian that includes the ef¬ 
fects of the Pauli term. Without either the mass 
term or the Pauli term, the theory is invariant 
under SU(2)l x SU(2)r chiral rotations, and its 
low momentum dynamics is described by the chi¬ 
ral Lagrangian 

C x = f -f Tr^S^S). (2) 

Here £ is an SU (2) matrix-valued field that trans¬ 
forms under the chiral group as £ —> LT,R\ with 
L are R being independent SU(2) rotations. Its 
vacuum expectation value, £o = (£), breaks the 
chiral symmetry down to an SU (2) subgroup. 
The fluctuations around £o correspond to the 
Goldstone-bosons, 

£ = £ 0 exp ji^TTaCTa/./V j . (3) 

The mass term explicitly breaks the underlying 
chiral symmetry, and a standard spurion analysis 
gives, to second order in to, the potential energy^ 

Vx = -jTr(£ + £t)+^{Tr(£ + £t)} 2 , (4) 

with Ci ~ to A 3 and C2 ~ to 2 A 2 . Since the 
Pauli term transforms under chiral rotations as 
does the mass term, and as it is proportional 

5 The symbol ~ means that we are not attempting to keep 
track of renormalization factors of order unity. 

6 Kinetic energy corrections to C x need not be considered 
since we are only interested in the vacuum state. 


to a single power of the lattice spacing, its ef¬ 
fects can be included (along with those of the 
mass term) by simply making the substitution 
m —> to + aA 2 , where A is required by dimen¬ 
sional analysis. Thus, we are led to the potential 
(f|) with coefficient^ 

ci ~ toA 3 + aA 5 , C 2 ~ m 2 A 2 + mat I 4 + a 2 A 6 .(5) 

We now study the alignment of the conden¬ 
sate, and the consequent pattern of symmetry 
breaking, as we reduce to at fixed (small) a. For 
a A 2 <C |to| <C A the discretization errors in ci and 
C 2 are a small correction, and symmetry break¬ 
ing is as in the continuum. In particular, since 
C 2 /C 1 ~ |to|/A <C 1, only the first term in the 
potential is important, and the condensate aligns 
as £0 = ±1. Flavor and parity are unbroken. 

This situation persists even when to becomes 
of size ~ aA 2 . Now discretization errors are im¬ 
portant, but C 2 /C 1 ~ aA is still small, and so 
the condensate is aligned as before. The effect 
of the discretization errors is to shift the mass at 
which ci vanishes from m = 0 to m' = 0, with 
m! = to — aA 2 (recall that all constants are be¬ 
ing set to unity). In terms of this new mass, the 
coefficients become 

ci ~ to'A 3 , C 2 ~ to' 2 A 2 + to'oA 4 + a 2 A 6 . (6) 

Finally, when the shifted masses are of 0{a 2 ), 
i.e. am' ~ (aA) 3 , the coefficients simplify to 

ci ~ m'A 3 , C 2 ~ a 2 A 6 . (7) 

The crucial point is that ci ~ C 2 , and, as we 
show below, competition between the two terms 
can lead to spontaneous flavor and parity break¬ 
ing, i.e. the appearance of an Aoki phascfj. 
It follows that the width of the Aoki phase is 
Atoo ~ a Am' ~ (aA) 3 . 

Let us now examine the potential (® in some 
detail to determine the condensate (£ 0 ). Writing 
£ = A + iB • cr with A 2 + B 2 = 1, the potential 
becomes 

V x = -ciA + c 2 A 2 , (8) 

7 Since we are keeping terms of OLo?) in C 2 , we must also 
include such terms in C e ff, Eq. (Dj). This has been done 
in [pi and does not change (p|). 

8 When the first two terms oft he potential are comparable, 
the cubic C 3 -term is down by a power of a A, except for a 
small region of order (aA ) 4 in which c\ ~ C 2 ~ C 3 . 




The parameter A is constrained to lie between — 1 
and +1 inclusive. Under flavor SU( 2)y rotations, 
A is invariant while B rotates by an orthogonal 
transformation. Hence, when the vacuum state 
Eo = Aq + *Bo • er develops a non-zero Bq, the 
flavor symmetry breaks spontaneously to a 17(1) 
subgroup defined by exp |z0Bo ■ er|. A non-zero 
value of Bo can occur only when Aq is strictly 
less than one. 

The sign of C 2 , which is not predicted by our 
analysis, determines two qualitatively different 
behaviors. We discuss c 2 > 0 first, since this gives 
rise to an Aoki phase. The potential is then a con¬ 
cave upward parabola with absolute minimum at 
A m = ci/2c 2 = e. If this minimum lies outside 
the range —1 to 1, then the vacuum is at Aq = ±1 
and flavor symmetry is not spontaneously bro¬ 
ken. But if \A m \ < 1, then the vacuum is at 
Aq = e, and flavor symmetry is spontaneously 
broken. The system thus finds itself in the Aoki 
phase with two massless Goldstone pions. The 
mass spectrum is 

m\ = = 0 , = 1 — e 2 for |e| < 1 (9) 

2c 2 

^ = |e|-lfor|e|>l, (10) 

which is illustrated in Fig. |2j. Recalling from (0) 
that e ~ m'/(a 2 A 3 ), we see explicitly that the 
Aoki phase has a width Amo ~ a Am' ~ (aA) 3 . 

The c 2 < 0 case gives a concave downward 
parabola with an absolute maximum at A m = e, 
and therefore the vacuum is always at Aq = ±1 
for any value of e. Flavor symmetry is not spon¬ 
taneously broken, and there are three degenerate 
pions with non-zero masses to 2 / 2 /2|c 2 | = 1 + |e|. 

In summary, in the continuum, when the quark 
mass passes through zero, the condensate jumps 
discontinuously from Eo = +1 to Eo = —1, with 
the pions becoming massless at the transition 
point. On the lattice, the explicit chiral symme¬ 
try breaking due to the Wilson term alters this 
transition. If c 2 > 0, the condensate rotates con¬ 
tinuously from +1 to —1, breaking flavor and par¬ 
ity as it does so. If c 2 < 0, the condensate jumps 
discontinuously as in the continuum, but the pi¬ 
ons never become massless. This shows how the 
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Figure 2. Pion masses as a function of e = Ci/2c 2 
for c 2 > 0. The curves are labeled by the flavor 
of the corresponding pion. 


Aoki phase is intimately related to chiral sym¬ 
metry breaking. Our analysis cannot determine 
the sign of c 2 , and indeed c 2 can have different 
signs for different choices of discretization of the 
Dirac operator. Numerical evidence to date is 
consistent with c 2 > 0 for both unimproved and 
improved Wilson fermions. 

We close with a few comments [||. First, one 
can relate the appearance of the Aoki phase to a 
non-vanishing density of zero modes of the her- 
mitian Wilson-Dirac operator. We stress, how¬ 
ever, that our analysis applies only to the un¬ 
quenched theory, and has no direct implications 
for the spectrum of the Wilson-Dirac operator in 
the quenched theory. Second, our analysis ap¬ 
plies as well to 0(a) improved Wilson fermions. 
And, finally, we stress that the phase structure 
will be different for Nf = 3, for which there is no 
m —> —m symmetry in the continuum. 
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